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We explored the anomalously low moment in MnSi found by experiment (~ 0.4 hb/Mti) vs the 
moment predicted by density functional theory (~ 1.0 /xb/Mu). With the addition of a Hubbard-U 
correction, we found several solutions with lower moments. These lower moment solutions show an 
unusual magnetic order and a magnetic quadrupole moment. The behavior of the moment under 
pressure does not follow the experimental trend. 
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I. INTRODUCTION 

A number of transition metal binary compounds can 
crystallize in the B20 s tructure and exhibit a wide variety 
of behaviors .1^^^'^'^'^ In particular, MnSi is a magnetic 
metal with an ordered moment n = 0.4 /zb/Mu, and 
the order manifests as a long period (~ 180 A) spin he- 
lix along the [111] direction due to the lack of inversion 
symmetry."- Above the relatively small Tc of 29 K, MnSi 
exhibits Curie- Weiss like behavior with a moment of 
fJ'CW = 2.2 /ie/Mn which indicates that the magnetism is 
due to itinerant electrons.f^ More recent attention has fo- 
cused on its metamagnetic transition with the application 
of ~ 1.46 GP a, above which MnSi exhibits non-Fermi liq- 
uid behavior .I^EESin] Similarly, FeGe is an itinerant heli- 
magnet with an ordered moment /i = 1.0 /ie/Fe, a longer 
period (~ 700 A), and a high temperature phase where 
the helix lies along the [100] direction which changes to 
the [111] direction in a transition with a large tempera- 
ture hysteresis (AT ~ 30 K)}'^ Additionally, the magnet 
order is likewise destroyed with the application of 20 GPa 
of pressure, but the high pressure phase shows a residual 
magnetization that is attributed to disorder.^i^ In con- 
trast, FeSi is a small gap semiconductor t hat und ergoes 
an unusual metal to insulator transition.'^IIIIini There 
is an ongoin g deba te on whether its prope rties a re due 
to Kondo-likd^l^E^ or band-like interactions,!^^!!!] among 
others'^^ 

Previous DFT calculations have had mixed results pre- 
dicting the behaviors of the B20 compounds. FeSi is 
correctly predicted to be a nonmagnetic semiconduc- 
tor with a gap that is comparable to experiment.!^ 
Similarly, FeGe is predicted to have a FM moment 
of 1.0 ^e/Fe at the experimental volume .I^^HHMls] 
However, MnSi has been predicted to have a similar 
moment between ~ 0.9 /xp/Mn and 1.2 yUe/Mn when 
calculated using the experimental latti ce const ant us- 
ing either an LMTO or LAPW method. ElasHlllll the 
LDA minimum volume, MnSi is predicted to be either 
paramagnetic!^ or ferro magn etic with a moment in the 
range 0.75 to 0.79 /Mn.l^ Calculations for both MnSi 
and FeGe have shown their metamagnetic transitions un- 
der pressure, but the predicted moment collapse requires 
at least twice the experimental pressure when starting 



from the experimental volumes The noted similar- 
ities between the predicted density of states of the B20 
monosilicides'^*' suggests that a fixed spin moment cal- 
culation may yield useful results,'2Sl but this solution is 
unsatisfactory as it adds an additional chemical poten- 
tial to the calculation while not addressing the underly- 
ing cause of the discrepancy. However, the narrowness 
of the Mn 3d - Si 3p bands compared to the Fe 3d - Ge 
4p bands and the non-negligible electronic interactions 
suggest that additional correlations are needed. '^^ 

In this paper, we examine the effect of adding a 
Hubbard-C/ interaction to our DFT calculations in the 
form of an LDA-I-U functional. In general, the Hubbard- 
U is a fixed property of the system. But, LDA-I-U lacks 
dynamic screening, so the calculated value usually dif- 
fers from the measured value.!^ Considering this and 
the inherent variability in DFT exhibited by the prior 
MnSi calculations, we chose to vary U over a range in- 
cluding the experimental value'^^ and a value calculated 
for MnCl^This allowed us to fit the experimental data 
better, and for values of U similar to the MnO value, 
we found a ground state with a moment comparable to 
experiment.!^ 

Depending on the value of U chosen, our calculations 
show two different ground states connected by a com- 
plex transition region. The low-U ground states retain a 
moment of 1 /ie/Mn while showing decreased magneto- 
volume coupling and becoming half-metallic as U is in- 
creased. However, the high-[/ ground states are very dif- 
ferent exhibiting a marked reduction in moment, a signif- 
icant rearrangement of the spin density, and a non-zero 
magnetic quadrupole moment around the manganese. 

Secti on [ll| details the specific DFT method we used. 
Section |III| explores the nature of the DFT ground state 
as a function of U and lattice constant. The appendix 
details how we calculated the magnetic quadrupole mo- 
ments. 



II. METHOD 

For our calculations we used the full potential, LAPW 
method within the local spin density approximation as 
implemented in Wien2k.'^'^ Each of our calculations were 
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performed on a 12 x 12 x 12 grid with 432 k-points in the 
irreducible Brillouin zone, a Fourier expansion of the po- 
tential up to a maximum reciprocal lattice magnitude of 
Gniax — 14, and muffin tin radii of Rmu — 2.197 a.u. and 
i?Si = 1-907 a.u.. The results did not vary with the size 
of the muffin tins up to a maximum of i?Mn = 2.3 a.u. 
and i?si = 2.0 a.u.. We set the core/valence separation 
energy to include the Si 3s and 3p and Mn 3p, 3d, and 4s 
electrons in valence. Also, we added local p and d orbitals 
to the Mn to eliminate ghost bands and improve conver- 
gence. The convergence criteria that we used was that 
the integrated absolute difference between the electron 
densities on subsequent iterations was less than 0.001 e. 
Additionally, to ensure that our calculations were well 
converged with respect to energy, we varied the number of 
plane waves from 269 to 1240 by varying the -Rmt *^max 
value between 5 and 9. To this, we added a Hubbard-C/ 
correction using the Around Mean Field (AMF) method 
to the 3d electrons on the manganese.*^ 

In the LDA-I-U AMF method, the following is added 
to the LDA energy functional for each atom and electron 
shell for which additional correlations are desired, 

m,m/,(T 

(1) 

where nma is the number of spin-cr in the m subband 
and flu is the average number of spin-cr present in the 
entire shell. The Coulomb, Ummi, and exchange, Jmmh 
matrices are defined in terms of Slater integrals, , and 
Clebsch-Gordan coefficients. These are connected to the 
Coulomb repulsion, U , and exchange, J, energies by set- 
ting U = F°,J= {F^+F-^)/U, andF^ = O.625F'^0We 
varied U over the range 0.0 Ry < U < 0.7 Ry, iterating to 
convergence at each value. Like U, The exchange energy, 
J, is a fixed parameter of the system, and from optical 
spectra of Mn embedded in a metallic matrix, it is esti- 
mated to be ~ 0.06 Ry."^^ Constrained LDA calculations 
for MnO give a similar result.!^ Our calculations show 
a small variation in the location of the transition region 
and the moment changes from /i ~ 0.245 to 0.173 /j-b/Mu 
as J is varied from 0.04 to 0.08 Ry at J7 = 0.60 Ry, and 
the results discussed in section IIIII use the atomic value 
of J = 0.06 Ry. 

Additionally, we added a spin-orbit term to give 
the system a preferred spin direction, stabilizing the 
Hubbard-U calculation. Initially, we chose [112] and [101] 
as our magnetization directions. However, they gave 
similar results to the magnetization along [100] direc- 
tion, which leaves the manganese ions equivalent, and 
is more computationally efficient. The results discussed 
within are for the magnetization along the [100] direc- 
tion. Finally, we simulated hydrostatic pressure by vary- 
ing the lattice constant between the experimental value 



of 4.558 A and 4.338 A, equivalent to a maximum pres- 
sure of nearly 20 GPa.™ 



III. RESULTS 
A. The Moment as a Function of U and Pressure 
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FIG. 1: (color online) The total moment per manganese (row 
1) and the percentage of the total moment in the interstitial 
regions (row 2) as a function of U at the both the experimental 
lattice constant (4.558 A, column a) and at high pressure 
(4.398 A, column b). The ground (excited) state solutions are 
represented by a filled red circle (blue square), and partially 
converged solutions are only the outlines. The three solution 
regions are labeled, and region II is highlighted. The arrows 
are guides to the eye indicating how U was changed along a 
particular branch. There are several solutions circled in (la) 
at the U values of 0.0, 0.4, 0.5, and 0.7 Ry, and they are the 
ones used below in the subsequent table and figures. 

In Fig. [l] we show the total magnetic moment and 
the moment found in the interstitial regions of MnSi 
as a function of U for two different lattice constants: 
the experimental lattice constant at atmospheric pressure 
(4.558 A) and at high pressure (4.398 A) corresponding 
to about 14.8 GPa. At a number of values of U, we found 
multiple solutions. We determined the ground state so- 
lutions by extrapolating from the total energy's variation 
with i?MT * -f'^max to determine a best estimate total en- 
ergy for each solution. In a couple of cases, the ground 
state could not be determined within the error of the fit, 
so both solutions were listed as the ground state. Ad- 
ditionally, the partially converged solutions are often in 
regions of function space where large rearrangements of 
charge and spin do not change the energy much, causing 
these solutions to require a large number of iterations to 
converge. 

We have divided our solutions into three regions: a 
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FIG. 2: (color online) The total moment as a function of 
lattice constant for U equal to (a) 0.0 Ry, (b) 0.4 Ry along the 
upper branch, and both (c) 0.6 Ry and (d) 0.7 Ry along the 
lower branch. The upper axis is the pressure as determined 
by using the experimental bulk modulus.'^ The gray vertical 
line is at the experimental metamagnetic transition point of 
1.46 GPa. 



high- moment, low-U region (region I); a low-moment, 
high-[/ region (III); and a complex transition between 
the two (II). At the experimental lattice constant, we 
found a single ground state which maintains a moment 
of 1 /zb/Mu in region I, and this moment is maintained 
into the transition region. In general, the moments of the 
solutions in region I decrease with pressure, as shown in 
Fig. |2] (a) and (b). But, this variation shrinks as U is 
increased, until U « 0.34 Ry where the moment does 
not change at the pressures we have looked at. Addition- 



ally, we hav e been able to reproduce the work of Yamada 
id Terao as below a lattice constant of ~ 4.336 A a 
non-magnetic solution becomes the ground state.^^^ 

Region I is also characterized by a highly localized spin 
density around the manganese atoms. Most of the excess 
spin density within the lattice is found within a distance 
of less than half of the muffin tin radius from the Mn. 
Additionally, as U is increased the percentage of the to- 
tal moment contributed by the interstitial regions grows 
from ~0.1%to~'5%at the region boundary. 

Following the 1 fiB /Mn moment solutions, we note that 
they quickly become metastable once region II has been 
entered, and the system switches to another solution with 
a lower moment as U is increased. At the experimental 
lattice constant, the system undergoes a transition to a 
state with a slightly reduced moment. This state does 
not seem to be accessible at any other lattice constant, 
and we consider it spurious. At higher pressure, the sys- 
tem undergoes a transition to a middle solution exhibit- 
ing hysteresis in U, and the moments of this branch are 
highly U dependent. With further increase in U, the 
system changes again to a lower moment solution that 
becomes the only accessible state in the third region. As 



U{Ry) 


LDA Min (A) % Difference 


0.00 (nonmagnetic) 


4.428 


-2.85 


0.00 


4.438 


-2.63 


0.40 


4.475 


-1.82 


0.60 


4.675 


2.56 


0.70 


4.619 


1.34 



TABLE I: The LDA minimum lattice constants were calcu- 
lated and compared to the experimental lattice constant for 
the solutions shown in Fig. [2] and the nonmagnetic solution. 



pressure is increased, the width of the region increases 
and shifts to the right. Also in this region, the contri- 
bution from the interstitial moment increases at a faster 
rate, approaching 16 % as region III is entered. 

The third region is characterized by a single low mo- 
ment state, with the exception of the spurious solution 
discussed above. The states with slightly higher moments 
have converged within our tolerance level, and should 
be considered the same solution as they would have the 
same moment if a more stringent criteria were imposed. 
The moments of the solutions in this region increase with 
pressure, as shown in Fig.[2](c) and (d). This is expected 
as the localizing effects of the Hubbard interaction are 
inversely proportional to the bandwidth, which is related 
to the interatomic spacing.^^^ So, as the pressure is in- 
creased, the system is pushed back towards the low-U 
ground states. Additionally, the interstitial contribution 
to the total moment remains high. But, at the experi- 
mental lattice constant it peaks at 18 % at ?7 « 0.6 Ry, 
and decreases to ~ 10 % at C/ = 0.7 Ry. 

Additionally, we examined the behavior of the LDA 
minimum lattice constant with U , as outlined in Tbl. [T] 
Curiously, it behaves differently for the low-U states and 
the high-?/ states. For the low-U states, the LDA min- 
imum is smaller than the experimental lattice constant 
and increase with U. This is expected behavior as LDA 
alone underestimates bond lengths and the additional re- 
pulsive interaction would force the atoms farther apart. 
However, the high-J7 LDA minimum lattice constants 
are larger than the experimental value and they decrease 
with increasing U. We believe this is caused by the de- 
creasing spatial overlap between the majority and mi- 
nority spin densities as U is increased, which reduces the 
pressure caused by the on site interaction between the 
majority and minority spins. 



B. Density of States 

In order to understand the nature of the transition 
from the low-U states to the high-JJ states, we looked 
at the density of states (DOS), as shown in Fig. |3] The 
most notable feature of our calculated DOS is the gap in 
the majority and minority spins. At U = Ry, there is 
a small gap in the minority spin, which disappears at the 
higher values of U. However, the majority gap increases 
in size and shifts towards Ei? as U increases. The high- 
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FIG. 3: The total (black line) and manganese d-shell (grey 
line) spin resolved density of states for all six of the circled 
solutions in Fig. [ija) are shown. The minority spin DOS is 
displayed along the negative vertical access. In (a), there are 
arrows pointing out the majority and minority gaps. 

moment solutions become half-metallic a.t U — 0.3 Ry, 
and they remain so until the system transitions to an- 
other state. For the low moment solutions, the gap re- 
mains large (> 1 eV) and has jumped above the Fermi 
energy. The second notable feature is the dramatic rear- 
rangement the DOS undergoes as U is varied. 

For a single s-band, the addition of a Hubbard- [/ splits 
the up and down states and forces them away from the 
Fermi energy.ES This essential feature is visible in the 
MnSi DOS as the d-shell states are forced away from each 
other as U is changed from 0.0 to 0.7 Ry. More precisely, 
at low values of U the manganese d-shell states make up 
nearly all of states available near the Fermi energy. But, 
as [/ is increased, the d-shell electrons are forced away 
from Fip flattening the peaks near Fip. For U = 0.5 Ry 
the Mn d states contribute about 70 % of the total weight 
at Ei?, and for U = 0.7 Ry this is reduced to less than 
50 %. This contributes to the dramatic rearrangement 
of the states that occurs with the increase in U, and the 
DOS of the interstitial states grows to fill in the majority 
of the states available at Ej? for high values of U. 

C. Angular Spin Distribution 

By integrating the DOS up to Ep, we see that the 
decrease in the moment is caused by a small transfer of 



(a) 



State 


0.0 (Ry) 0.4 (Ry) 


0.5 (Ry) 0.7 (Ry) 


^x-^ —y'^-\-xy 
^xz-\-yz 


0.21 0.22 
0.40 0.41 
0.39 0.37 


0.18 0.22 
0.39 0.41 
0.43 0.38 


Total 


3.05 2.99 


2.62 2.60 


(b) 


State 


0.0 (Ry) 0.4 (Ry) 


0.5 (Ry) 0.7 (Ry) 


^2z'^—x'^-y'^ 
^x^ —y'^-\-xy 
^xz-\-yz 


0.18 0.11 
0.39 0.37 
0.43 0.52 


0.20 0.13 
0.40 0.38 
0.40 0.49 


Total 


2.03 2.08 


2.44 2.46 



TABLE II: The relative occupancies of the Mn d-subshells 
and the total occupancy of the Mn d-shells are shown for the 
majority (a) and minority (b) spins for several values of U . 
The tables are divided into the low-f/ and high-C/ branches. 

electrons from the majority to minority spin d-bands, as 
shown in Tbl. HT] This also reveals that the transfered 
electrons are not evenly distributed among the minority 
spin bands. While the rearrangement is not large, it does 
cause a visible shift in the alignment of the spin density 
around the Mn atoms, as shown in Figs. [4]and[5j 

The most notable feature of the spin densities is the 
minority spin trilobal structures that are centered around 
the [111] direction. Within in any branch of the solutions, 
the trilobal structure is distorted slightly as U is varied. 
However, the degree of spin polarization increases with 
increasing U in the high moment branch (Figs. |4]^a), (b) 
and[5|. In the case of the low-moment solutions, some 
regions have saturated and the boundary regions between 
the spin polarizations decrease in size. Additionally, the 
center of the trilobal structures changes spin polarization 
as U is increased in the low-moment branch. Between the 
two solution branches, the trilobal structure is rotated 
about the [111] direction which tends to align the lobes 
with six of the thirteen nearest neighbors. 

This arrangement of spins is not unique to the hav- 
ing the magnetization aligned along the [100] direction. 
With the magnetization aligned along the [101] direction 
at U = 0.5 Ry, the arrangement of the spin density is 
nearly the same as that shown in Fig. ^c). Since, the 
structure does not seem dependent upon the magneti- 
zation direction, nor do we find that it is affected by a 
10 T field, the anisotropy may be detectible as a small 
magnetic quadrupole moment around the Mn atoms. 

D. Magnetic Quadrupole Moments 

From Ref.|39l the magnetic quadrupole due to a non- 
uniform spin density can be calculated from the traceless 
tensor operator 

Mtj = TiSj + rySi - 25ijrkSk (2) 
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FIG. 4: These are plots of the spin densities around the Mn 
atom radially integrated out to the muffin tin boundary for 
the circled solutions along the upper and lower branches in 
Fig. [TJa). The azimuthal angle, 6, is plotted on the vertical 
axis and the polar angle, (p, is plotted on the horizontal axis. 
The directions to the nearest neighbor Mn (Si) atoms are 
plotted as black circles (gray squares). 




4.392 and 4.558 A as a function of U. 

For the low-U ground states, the magnetic quadrupole 
moments are negligible, and are further suppressed by 
the increase in pressure. In the transition region, the 
quadrupole moments increase quickly, and show hystere- 
sis with U that was seen in the dipolc moments. The 
spurious solution at the experimental lattice constant, 
increases to the value obtained by the middle solution at 
U = 0.46 Ry, but it does not increase beyond that point. 
The middle solution quadrupole moments continue to in- 
crease as U is increased. The quadrupole moments are 
nearly their largest value of 0.15 /ZBa.u. as the high-C7 re- 
gion is entered. At the experimental lattice constant, the 
largest magnetic quadrupole is found at ^ 0.58 Ry, and 
it decreases with further increase in U. Since an increase 
in pressure shifts the transition region to the right, the 
quadrupole moments continue to increase at least until 
U = 0.7 Ry for a lattice constant of 4.392 A. Addition- 
ally, the quadrupole moments around the silicon atoms 
does not exceed 0.03 /iBa.u., and the quadrupole mo- 
ments in the nonmagnetic solution never exceeds a tenth 
of that value. 

While these moments are small, the principal direc- 
tions of the magnetic quadrupole are not aligned with 
the electric field gradient (EFG) at the Mn nucleus. 
The angle between the eigenvectors for largest magnetic 
quadrupole moment and the principal direction of the 
EFG is about 27° for U between 0.5 and 0.6 Ry, and over 
that range, the orientation of the magnetic quadrupole 
stays nearly constant. The principal direction of the 
EFG, also, remains roughly constant over that range, 
but the other directions are rotated ~ 24° about the 
principal direction. It should be noted, that the princi- 
pal directions of the magnetic quadrupole moment are 
mostly insensitive to pressure, and while the principal 
direction of the EFG is also insensitive to pressure, the 
other EFG directions rotate about 21° from their origi- 
nal positions. Since both Mossbauer spectroscopjEl and 
^^Mn NMlPSlare capable of detecting anisotropies in the 
hyperfine field, they should be able to detect the mag- 
netic quadrupole moment. 



(a) U = 0.00 Ry (b) U = 0.70 Ry 

FIG. 5: This is the same data as Fig. [4] (a) and (d) plotted 
on the surface of a sphere the size of the Mn muffin-tin. In- 
cluded as a visual reference are the nearest neighbor atoms 
with lighter (darker) shading for the Si (Mn) atoms. The 
bonding Si atom is labeled. 

where and Si are the ith components of the posi- 
tion and spin operators, respectively. In the appendix, 
we show this to be equivalent to the classical mag- 
netic quadrupole.'^ Fig. |6] shows the two independent 
principal components of the magnetic quadrupole mo- 
ments within the Mn muffin-tins for lattice constants 



IV. DISCUSSION 

Empirically, U was found to be 0.13 Ry, assuming that 
J = 0.06 Ry and neglecting the multipole contribution.!^ 
At this value of U, the moment remains 1 //s/Mn and it 
is less sensitive to pressure, so this value of U is clearly 
unsuitable. In an LDA calculation for MnO, a value of 
U = 0.51 Ry was found.l^The ground state at the com- 
puted value lies on the lower branch within the transi- 
tion region with a moment of ~ 0.24 /ib/Mu. While 
this solution has a moment that is comparable to ex- 
periment, the hyperfine field around the Mn nucleus is 
Hup — 13 kG when the orbital and dipolar contribu- 
tions are included. However, increasing U to 0.62 Ry 
gives = — 2.1 kG, which is in good agreement with 
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FIG. 6: The two independent magnetic quadrupole moments 
as a function of U for the lattice constants (a) 4.558 A and (b) 
4.392 A are shown. The larger (smaller) moment is labeled 
by circles (triangles) and squares (diamonds) and lies above 
(below) the x-axis. As in Fig. [l] the three solution regimes 
are highlighted and labeled and the arrows serve the same 
purpose as before. The quadrupole moments of the spurious 
solution are not shown. 



experiment.^ The hyperfine field at the Si nucleus, how- 
ever, is about half the measured —91.3 kOe and aligned 
parallel with the magnetic momentj^^ but the moment 
is only 0.19 /ie/Mn. As for the pressure dependence 
of these low moment states, we know that the experi- 
mental high pressure state exhibits a complex magnetic 
order.l^ Since the magnetic moment is well reproduced 
with LDA-I-U, this suggests that dynamic correlations are 
responsible for the observed pressure dependence. 

In summary, we have found a set of solutions with a 
moment comparable to experiment for values of U that 
are in line with prior work on other Mn compounds. 
These solutions exhibit a distinct magnetic anisotropy 
that manifests as a small magnetic quadrupole moment 
around the Mn, which should be experimentally de- 
tectable. Additionally, they correctly determine the hy- 
perfine field at the Mn atom. However, these solutions 
behave incorrectly under pressure and predict the wrong 
hyperfine field at the Si atom. 



Acknowledgments 

We thank John DiTusa and Richard Kurtz for useful 
conversations. This work is supported in part by NSF- 



OCI Grant No. 043812 and a Faculty Research Grant 
from the LSU Office of Sponsored Research. 

APPENDIX: CALCULATING MAGNETIC 
QUADRUPOLE MOMENTS 



According to Ref.|4Di the multipole moments of the 
magnetization density, M(x), can be calculated by the 
formulations for the moments of the charge density by 
substituting pn — V • M for the charge density. The 
cartesian form of the magnetic quadrupole is then 



s. 



(M) ^ 
ij 



{SxiXj - 6ijXixi)dkMk d^x, (A.l) 



where the summation convention has been employed and 
di — d/dxi. To connect this formula to the operator in 
Eqn. ([2]), we must integrate by parts. By first noting that 
diXj — Sij, the first term in the integrand becomes 



XiXjdkMk = dkixiXjMk) - (xiMj + XjMi) (A. 2) 



and the second term is 



SijXixidkMk = 5ij[dk{xixiMk) - 2xiMi]. (A.3) 



Substituting these into Eqn. (A.l I, we get 



Q-^^ \ I dk{xiXjMk - SijXiXiMk)d^x 

+ ^ I [3{xiMj + XjM,) - 26^jXkMk] d^x (A.4) 



The first term is the surface term, and in the case of a 
localized density it can be taken to be zero.l^ Since we 
are integrating over the Mn muffin tins, the density is 
not localized and there are spurious monopole terms if 
the surface terms are ignored. The second term is the 
expectation value of the magnetic quadrupole operator 
m Eqn. This can be seen by noting that the spin 
density, = P] — P[, is related to the ground state 
spinor, ^ , by ■ (7^ -^ where is the Pauli matrix. In 
our calculations, we use expectation value of M directly, 
as it contains the surface term. 
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